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DOPPLER  LIDAR  FEASIBILITY  STUDY 
FOR  THE  BALLISTIC  WINDS  INITIATIVE 


1.  INTRODUCTION 


1.1.  Background  and  Purpose  of  This  Study 

The  goal  of  the  Ballistic  Winds  (BW)  initiative  is  to  improve  the  accuracy  of  bomb  drops 
from  high  altitude  aircraft.  It  has  been  shown  that  knowing  the  vertical  wind  profile  can 
markedly  improve  the  accuracy  of  ballistic  weapons.  As  part  of  the  BW  program,  the  high 
altitude  wind  field  will  be  measured  with  an  onboard  coherent  lidar  system.  The  measurements 
will  then  be  used  by  the  onboard  targeting  system  to  make  more  refined  real  time  wind 
corrections. 

The  current  study  investigates  the  feasibility  of  obtaining  the  wind  field  profile  for  real 
time  corrections  during  targeting  nms.  This  study  was  primarily  concerned  with  obtaining  a 
methodology  for  predicting  the  performance  of  coherent  lidar  systems  which  could  be  used  in 
both  the  design  and  field  test  measurement  phases  of  the  Ballistic  Winds  Initiative  Program. 

Some  Ballistic  Winds  field  tests  are  planned  to  utilize  a  B-52  aircraft  fitted  with  a 
Lightwave  Laser  Inc.  tunable  Tm:YAG  solid  state  laser  Doppler  lidar  system  at  a  fixed  off  nadir 
look  angle  and  locked  in  a  fixed  azimuthal  direction.  The  proposed  aircraft  will  fly  in  a 
rectangular  racetrack  flight  path  at  a  test  range  with  the  bomb  drop  made  in  the  north-south 
direction  as  shown  in  Figure  1.  The  tests  will  occur  during  the  late  winter  to  early  spring  time- 
frame.  For  the  White  Sands  Missile  Range  (WSMR)  test  site,  this  is  when  there  are  prevailing 
westerly  winds.  The  initial  bomb  drop  will  use  the  standard  correction  procedure  using  only  the 
aircraft  altitude  wind  speed  as  input.  However,  the  lidar  will  take  measurements  for  later 
analysis. 
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Figure  1.  Nominal  flight  path  for  Ballistic  Winds  lidar  tests. 

1.2.  Organization  of  the  Report 

Section  2  describes  state  of  the  atmosphere  at  a  mid-latitude  site  in  spring  including 
expected  wind  conditions  and  the  atmospheric  attenuation  and  backscatter  due  to  the  molecular 
and  aerosol  constituents.  The  atmospheric  models  including  the  aerosols,  are  taken  from 
LOWTRAN?^  and  are  only  “nominal”,  average  conditions.  They  do  not  give  a  true 
representation  of  actual  site  conditions  and  atmospheric  measurements  should  be  made  at  any 
site  planned  for  BW  tests  including  WSMR.  Section  3  describes  the  lidar  simulation  parameters 
used  in  the  feasibility  analysis  and  the  signal-to-noise  (SNR)  performance  relations  along  with 
the  methodology  used  in  the  prediction  analyses.  Section  4  describes  the  lidar  simulation  SNR 
performance  predictions.  This  section  includes  a  discussion  of  the  utility  of  pulse  integration  to 
improve  the  SNR,  and  range  and  velocity  error  considerations.  Section  5  gives  a  brief  discussion 
of  the  wind  effects  on  the  aircraft  and  velocity  measurements.  Summary  and  conclusions  are 
given  in  Section  6. 


^  Kneizys,  F.X.  Shettle,  E.P.,  Abreu,  L.W.,  Chetwynd,  J.H.,  Anderson,  G.P.,  Gallery,  W.O., 
Selby,  J.E.A.,  and  Clough,  S.A.  (1988)  Users  Guide  to  LOWTRAN7,  Air  Force  Geophysics 
Laboratory,  Hanscom  AFB,  MA,  AFGL-TR-88-0177,  ADA  206773. 
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2.  ATMOSPHERIC  EFFECTS 


When  evaluating  the  performance  of  a  lidar  system,  the  state  of  the  atmosphere  must  be 
included  in  the  simulation.  In  analyzing  the  signal  return  to  the  lidar  receiver,  all  the  atmospheric 
effects  that  degrade  and  influence  system  performance  must  be  taken  into  account.  Beam 
attenuation  comes  from  molecular  scattering  and  absorption,  and  aerosol  absorption  and 
scattering.  The  aerosol  backscatter,  which  gives  the  lidar  system  its  return  signal,  will  also  be 
discussed.  For  a  Doppler  system,  the  motion  of  the  aerosols  is  assumed  to  be  representative  of 
the  wind  velocity  and  direction. 

The  atmospheric  conditions  assumed  for  the  majority  of  the  analyses  done  for  this  study 
are  a  mid-latitude  spring/summer  model  with  low  humidity  (40%)  in  both  the  boundary  layer  and 
troposphere.  The  model  tropopause  boundary  can  vary  from  10  to  12  km,  depending  on  the 
location  of  the  jet  stream.  Background  stratospheric  conditions  from  LOWTRAN7  are  assumed. 
The  top  of  the  boundary  layer  is  at  2  km  with  a  23  km  surface  visibility  and  a  rural  aerosol 
model.  This  later  region  is  not  important  since  below  3  km  altitude,  a  bomb  is  supersonic  and 
crosswinds  do  not  affect  its  trajectory. 

Attention  must  be  paid  to  the  results  obtained  from  a  model  atmosphere  when  making 
system  predictions.  The  results  do  not  always  portray  the  true  conditions  at  a  test  site  and  the 
atmospheric  constituents’  loadings  can  vary  by  orders  of  magnitude.  These  variations  can 
greatly  affect  the  conclusions  drawn  from  predictions  and  care  should  always  be  exercised. 

2.1.  WSMR  Wind  Field 

WSMR  is  just  one  of  the  areas  being  considered  for  BW  field  tests.  It  has  limited  access 
and  the  atmosphere  has  been  measured  and  characterized.  The  aircraft  will  fly  at  the  edge  of  the 
jet  stream,  between  10  and  12  km  altitude.  At  these  altitudes,  the  jet  stream  winds  at  WSMR 
vary  between  80  and  120  kts  (~41  to  62  m/s).  The  WSMR  area  monthly  average  wind  profiles 
for  February  through  April  are  shown  in  Figure  2.  These  profiles  were  obtained  from  an  NCAR 
climatology  database^  for  the  1950-1964  time  period. 

2.2.  Atmospheric  Attenuation 

The  beam  attenuation  from  molecules  at  the  laser  wavelengths  of  interest  is  significant. 
The  Lightwave  Laser  system,  a  Tm:YAG  solid-state  laser,  is  tunable  from  2.0133  to  2.0138  pm. 
Using  FASCOD3P,^  Figure  3  shows  the  vertical  transmission  due  to  absorption  covering  the 
tunable  laser’s  wavelength  region.  The  propagation  path  is  from  12.5  to  3  km  for  mid-latitude 


^  Jenne,  R.L.  (1994)  Jenne’s  Northern  Hemisphere  Climatology  -  Monthly  (1950  -  1964)  Data 
Set  DS205.0,  NCAR  Research  Data  Archives,  Available  over  the  Internet 
(ftp;//ncardata.ucar.edu/datasets/ds205.0). 

^  Anderson,  G.P.,  Chetwynd,  J.H.  (1993)  “Fast  Atmospheric  Signature  Code  Version  3  - 
Preliminary”,  Presented  at  The  Joint  Phillips  -  Wright  Laboratories  Atmospheric  Propagation 
Workshop,  Wright-Patterson  AFB,  OH,  18-20  May,  Handout  #8. 
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Figure  3.  Vertical  transmission  from  12.5  km  to  3  km  altitude  due  to  molecular  absorption.  This 
is  for  propagation  within  the  2.0130  to  2.0140  pm  region. 
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Altitude  (km) 


spring  and  fall  conditions.  The  curve  also  shows  the  locations  of  the  major  molecular  absorbers 
(CO2,  H2O,  and  N2O).  The  minimum  wavelength  of  the  Lightwave  Laser  system  {i.e.,  2.0133 
jxm)  is  near  the  peak  of  the  transmission  curve  (2.0132  pm). 

In  this  study,  future  references  to  the  “peak”  line  refer  to  2.0132  pm,  while  the  “off  peak” 
line  refers  to  the  best  atmospheric  transmission  (2.0133  pm)  for  the  Lightwave  Laser  system. 

The  molecular  absorption  profiles  for  the  “peak”  and  “off  peak”  laser  lines  are  shown  in 
Figure  4  for  typical  mid-latitude  atmospheres.  These  values  were  obtained  fi'om  FASCOD3P. 
The  sharp  increase  in  molecular  absorption  below  3  km  is  due  to  the  waver  vapor  continuum. 

The  extinction  coefficient  profiles  for  the  “peak”  laser  line  are  shown  in  Figure  5.  The 
aerosol  component  is  based  on  the  standard  Geophysics  Laboratory  models'*  and  the  other 
atmospheric  conditions  defined  previously.  The  total  extinction  is  due  almost  entirely  to 
molecular  effects.  At  high  altitudes,  the  molecular  extinction  is  over  two  orders  of  magnitude 
greater  than  that  due  to  aerosols.  Molecular  attenuation  is  due  primarily  to  absorption  as  shown 
in  Figure  6. 


Figure  4.  Molecular  absorption  profiles  for  mid-latitude  atmospheric  models.  These  profiles  are 
for  “peak”  and  “off  peak”  laser  lines. 


Fenn,  R.W.,  Clough,  S.A.,  Gallery,  W.O.,  Goode,  R.E.,  Knei2ys,  F.X.,  Mill,  J.D.,  Rothman, 
L.S.,  Shettle,  E.P.,  and  Volz,  F.E.  (1985)  “Optical  and  Infrared  Properties  of  the  Atmosphere”, 
Chap  18  in  Handbook  of  Geophysics  and  the  Space  Environment,  A.S.  Jursa  (Scientific  Editor), 
Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA,  AFGL-TR-85-0315,  ADA  167000. 
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Figvire  5.  Extinction  coefficient  profiles  for  the  “peak”  laser  line. 


Figure  6.  Attenuation  coefficient  profiles  at  the  “peak”  laser  line. 
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2.3.  Backscatter  Coefficients 


The  profiles  of  the  aerosol  and  molecular  backscatter  coefficients  are  shown  in  Figure  7. 
These  profiles  are  valid  over  the  entire  wavelength  region  of  interest.  The  molecular  backscatter 
can  only  further  degrade  the  velocity  measurements.^  There  is  significant  contributions  to  the 
backscatter  coefficient  at  short  wavelengths  from  molecules.  This  effect  can  mask  any  base 
velocity  due  to  winds.  Because  the  molecules  are  traveling  randomly  with  a  mean  velocity  near 
the  speed  of  sound,  the  energy  returned  from  the  backscattered  molecules  is  randomly  spread 
over  a  wide  Doppler  spectrum  with  little  energy  within  each  velocity  bin.  Therefore,  the  Doppler 
shift  from  these  molecules  cannot  be  used  to  estimate  wind  velocity.  In  the  signal-to-noise 
calculations,  which  will  be  discussed  later,  molecular  extinction  due  to  both  absorption  and 
scattering  is  included,  but  the  backscatter  for  the  coherent  system  only  includes  the  aerosol 
backscatter  signal.  It  will  therefore  have  little  effect  on  the  system  signal-to-noise. 


Figure  7.  Profiles  of  the  molecular  and  aerosol  backscatter  coefficient. 


^  S.E.  Moody  (1987)  “Evaluation  of  Laser  Technologies  for  On- Aircraft  Wind  Shear  Detection,” 
SPIE  Vol.  783,  p.l24. 
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2.4.  Line-of-Sight  Transmission 


The  line-of-sight  (LOS)  transmission  losses  can  be  calculated  within  BACKSCAT  4.0^ 
and  LOWTRAN7,  but  simplified  transmission  calculations  were  made  to  quickly  get  a  sense  of 
the  variation  due  to  altitude,  season,  viewing  angle,  and  laser  line.  The  relations  between  the 
extinction  coefficient  and  the  vertical  distance  traversed  make  use  of  the  geometry  shown  in 
Figure  8. 


Figure  8.  Line-of-sight  transmission  geometry. 


The  vertical  transmission  (Xy)  is  defined  as 


Ty  =  exp 


—  Ja(5)<is' 
\  h,  j 


(1) 


where  a  is  the  extinction  coefficient  and  s  is  altitude.  Thus,  the  LOS  transmission  (x)  is. 


^  Longtin,  D.R.,  Cheifetz,  M.G.,  Jones,  J.R.,  and  Hummel,  J.R.  (1994)  BACKSCAT  Lidar 
Simulation  Version  4.0:  Technical  Documentation  and  Users  Guide,  Phillips  Laboratory, 
Directorate  of  Geophysics,  Hanscom  AFB,  MA,  PL-TR-94-2170,  ADA285851. 
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T  =  exp 

-  ^aiR)dR 

? 

\  R2 

(  *1 

«  exp 

-  Ja(5)  sec(p;yJ5  , 

V  h2 

j 

^  V 

(2) 


The  expressions  in  Eq.  (2)  are  approximate  because  a  plane-parallel  atmosphere  is  assumed. 

The  “peak”  laser  line  one-way  cumulative  molecular  transmission  at  30°  off  nadir  shown 
in  Figure  9  for  three  different  aircraft  cruising  altitudes.  By  3  km  altitude,  there  is  not  much 
difference  between  the  cumulative  transmissions,  less  than  2%.  The  off  nadir  vievdng  angle 
variation  for  the  “peak”  laser  line  with  aircraft  altitudes  of  12  km  and  10  km  is  shown  in 
Figure  11.  Figure  10  shows  the  effect  of  season  on  the  cumulative  molecular  transmission  for 
the  “peak”  and  “off  peak”  laser  lines.  Although  there  are  considerable  losses  of  up  to  44%,  the 
variations  shown  on  individual  figures  are  not  large. 


Figure  9.  Effect  of  aircraft  cruising  altitude  on  one-way  cumulative  molecular  transmission. 
These  transmission  profiles  are  for  the  “peak”  laser  line  only. 
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Figure  10.  Effect  of  off  nadir  viewing  angle  on  one-way  cumulative  molecular  transmission. 
These  are  for  the  “peak”  laser  line  from  aircraft  altitudes  of  (a)  12  km  and  (b)  10  km. 
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Figure  11.  Laser  line  one-way  cumulative  seasonal  molecular  transmission.  These  profiles  are 
for  “peak”  and  “off  peak”  laser  line  transmissions.  FW  =  fall/winter,  SS  =  spring/summer. 


3.  BALLISTIC  WINDS  LIDAR  SIMULATION  METHODOLOGY 


3.1.  Ballistic  Winds  Simulation  Parameters 

This  section  describes  the  system  and  simulation  parameters  used  in  the  performance 
study.  The  parameters  deal  with  the  laser  transmitter,  the  optics,  and  the  receiver/detector  values. 
The  B-52  aircraft  flies  between  10  and  12  km  altitude  above  the  ground  with  a  speed  of  220  to 
225  m/s.  The  lidar  is  fixed  at  an  off  nadir  viewing  angle  of  either  30°  or  most  likely  15°  due  to 
mounting  and  viewing  limitations  in  the  aircraft  fuselage. 

The  Lightwave  Laser  system  laser  transmitter  has  the  following  specifications  - 

1.  Tunable  solid-state  Tm:YAG 

2.  Wavelength:  X  =  2.0133  pm  (“off  peak”)  to  2.0138  pm,  “peak”  laser  line  of  2.0132 
pm  also  used 

3.  Pulse  energy:  Ep  =  3  mJ 

4.  Pulse  length:  tp  =  0.25  ps 

5.  Pulse  repetition  frequency  (PRF):  v  =  180  Hz 
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The  laser/aircraft  optics  has  the  following  parameters  - 

1 .  Aperture  diameter:  7  cm 

2.  Overall  lidar  optical  efficiency  (two-way,  transmit  &  receive):  =  10% 

3.  A/C  quartz  window  transmission  (two-way):  =  94% 

The  receiver/detector  parameters  are  - 

1.  APD  detector:  EPITAXX  lOOT  GR  2.2 

2.  Quantum  efficiency:  q  =  61% 

3.  Excess  noise  factor:  F=  1.2 

4.  Detector  system  bandwidth:  5  =  2  MHz. 

3.2.  Method  of  Evaluating  the  SNR  Performance 

The  SNR  performance  predictions  were  made  using  BACKSCAT  Version  4.0  in  the 
pulsed  coherent  Doppler  lidar  simulation  mode.  The  molecular  absorption  was  included  by 
independently  running  FASCOD3P,  creating  absorption  profiles,  and  then  inputting  them  into 
BACKSCAT.  The  molecular  backscatter  was  not  included  in  the  return  power,  although 
attenuation  due  to  molecular  scattering  was  included.  The  user  defined  detector  option  is  used 
for  the  EPITAXX  avalanche  photodiode  (APD).  The  outputs  from  BACKSCAT  used  in  this 
study  are: 

1 .  the  lidar  return  vs.  altitude  and  range, 

2.  the  single  pulse  SNR  (voltage)  vs.  altitude  and  range,  and 

3.  the  velocity  and  range  accuracy  vs.  altitude  and  range. 

BACKSCAT  runs  were  used  to  predict  the  single  pulse  signal-to-noise  performance  of  a 
coherent  Doppler  lidar  system  looking  towards  the  ground  from  the  aircraft.  Variations  in 
molecular  attenuation,  lidar  vievdng  angle,  A/C  altitude,  and  laser  line  propagation  were 
considered.  The  amount  of  pulse  summing  required  to  increase  the  SNRp  to  10  dB  was 
examined.  The  results  of  these  calculations  "will  be  discussed  in  Section  4.0. 

3.3.  Simulation  Assumptions 

In  BACKSCAT,  the  lidar  transmit  beam  is  assumed  to  be  collimated  and  there  is  no 
focusing.  The  transmit  and  receive  optics/beams  are  collinear  which  is  a  monostatic  system  with 
no  optics  offsets.  The  local  oscillator  power  is  assumed  to  be  large  enough  such  that  the  system 
is  shot  noise  limited  with  the  SNR  dominated  by  the  local  oscillator  noise  source. 

In  this  preliminary  study  there  are  no  turbulence  or  speckle  effects  included  which  can 
degrade  the  beam  and  reduce  the  return  signal.  The  local  oscillator  and  receive  beam  phase 
fronts  are  assumed  to  be  matched  so  that  there  are  no  phasing  losses. 
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3.4.  Equations  in  BACKSCAT  4.0 


The  voltage  signal-to-noise  (SNRy),  is  calculated  for  a  single  pulse  as 


SNR, 


'  nPs 
hvBF’ 


(3) 


where  Pj  is  the  return  power  into  the  receiver,  ti  is  the  detector’s  quantum  efficiency,  h  is  Planck 
constant,  v  is  the  laser  frequency  (=  c/X),  c  is  the  speed  of  light,  F  is  the  detector  excess  noise 
factor,  and  B  is  the  bandwidth  of  the  receiver  electronics  assumed  to  be  a  matched  filter. 

The  return  signal  power  is  defined  as 


=  E^TTXpxix^ 


.2 

=  1.63x10’'  " 


km 


(4) 


where  p  is  the  aerosol  backscatter  coefficient  in  m'^  sr'*,  is  the  receiver  collector  area,  x^  is  the 
one-way  atmospheric  transmission,  Xj  is  the  transmitter  optical  transmission,  Xr  is  the  receiver 
optical  transmission,  x^  is  the  one-way  transmission  through  the  quartz  window  on  the  aircraft, 
and  R  is  the  range  from  the  aircraft  to  a  given  altitude  at  the  specific  viewing  angles  of  the  lidar. 

Substituting  in  for  the  lidar  system  parameters,  the  SNRy  becomes. 


SNR,  =  1.61  X  10^  7^, 

=  2.04xl0^^f^. 

In  the  SNR  relation,  it  is  assumed  that  a  matched  filter  is  used  in  the  detector  electronics,  i.e., 


(5) 


(6) 


13 


4.  SNR  PERFORMANCE  PREDICTIONS 


4.1.  Single  Pulse  Operation 

BACKSCAT  runs  were  made  to  predict  the  various  effects  that  can  influence  the  single 
pulse  operation  of  the  Ballistic  Winds  lidar.  These  results  are  shown  in  Figures  12  through  15 
where  the  molecular  effects,  the  viewing  angle,  the  aircraft  altitude,  and  the  laser  lines  were 
individually  varied.  In  all  of  these  calculations,  molecular  backscatter  is  not  included  in  the 
signal  return,  except  for  the  investigation  of  the  various  molecular  effects  (Figure  12). 

The  trends  in  each  case  are  similar,  the  voltage  SNR  is  relatively  large  near  the  aircraft, 
enough  for  single  pulse  detection  of  the  return.  The  signal  decreases  due  to  the  increasing  range 
with  insufficient  backscatter  to  compensate  for  the  range  losses  until  near  the  boundary  layer. 
Although  the  aerosol  return  is  largest  in  the  boundary  layer,  this  region  does  not  affect  the 
ballistic  trajectory  and  the  returns  from  this  region  would  not  be  processed  in  any  operational 
system.  When  all  the  various  scenarios  analyzed  are  compared,  there  is  not  a  great  variation  in 
the  resulteint  SNR. 

This  indistinguishable  results  in  Figure  15  are  especially  important  for  the  laser  line  SNR 
predictions.  That  is,  the  Lightwave  Laser  system  is  not  able  to  tune  itself  to  the  peak  of  the 
molecular  transmission  line,  but  the  effect  of  being  “off  peak”  is  unimportant. 


Altitude  (km) 

Figure  12.  Molecular  attenuation  variation  effects  on  BACKSCAT  single  pulse  SNR. 
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Single  Pulse  SNR,  Single  Pulse  SNR 


Altitude  (km) 


Figure  15.  Laser  line  variation  effects  on  BACKSCAT  single  pulse  SNR. 


4.2.  Pulse  Integration  Performance 

Pulse  integration  can  increase  the  signal-to-noise  by  summing  single  pulses  such  that  the 
return  signal  power  is  increased  at  a  faster  rate  than  that  of  the  noise  power.  Ideal  pulse  also 
called  coherent  integration,  increases  the  resultant  SNR  by  the  number  of  pulses  summed.  In 
incoherent,  or  post  detection  integration,  the  signal-to-noise  is  increased  by  the  square  root  of  the 
number  of  pulses  summed.  With  most  pulse  integration  processes,  the  integrated  signal-to-noise 
(SNRj)  falls  between  coherent  pulse  integration  with  a  gain  (y)  of  one  and  post  detection 
integration  with  a  gain  of  0.5,  i.e., 

SNRi  =  iV^  SNRi  (0.5  <  y  <  1).  (7) 

where  A  is  the  number  of  pulses  summed  and  SNRj  is  the  single  pulse  signal-to-noise  ratio. 

There  are  many  factor  which  affect  the  integration  efficiency.  These  include  the  number 
of  pulses  summed,  the  receiver  electronics,  the  target  characteristics,  and  the  pulse  timing,  shape, 
etc.  As  a  worst  case,  it  is  assumed  that  the  ballistic  winds  pulse  summing  process  is  incoherent 
(y  =  0.5). 
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When  summing  the  lidar  return  pulses,  it  is  assumed  that  the  wind  field  remains  uniform 
and  constant  over  the  number  of  pulses  required  to  reach  the  desired  integrated  SNR.  The  time 
limitations  to  this  assumption  will  have  to  be  determined  during  the  field  tests. 

In  the  following  analysis,  it  is  assumed  that  10  dB  power  SNR  is  required  for  wind  field 
detection  in  each  range  bin.  This  SNR  (power)  is  equivalent  to  a  SNR  (voltage)  of  5  dB 
(equivalent  to  a  ratio  of  3.2).  Figures  16  through  19  show  the  required  number  of  integrated 
pulses  necessary  to  reach  this  total  return  power.  The  figures  are  in  the  same  order  as  the  single 
pulse  SNR  figures  shown  previously.  For  an  aircraft  flying  at  10  km,  the  maximum  number  of 
pulses  that  must  be  summed  to  maintain  10  dB  signal-to-noise  (power)  is  64.  This  occurs  in 
approximately  the  5.5  km  altitude  region. 

Figure  20  shows  the  required  number  of  pulses  that  must  be  summed  to  maintain  a 
desired  integrated  SNR  while  the  various  molecular  effects  are  accounted  for.  This  is  for  an 
aircraft  cruising  at  10  km  altitude  using  the  “peak”  laser  line. 


Altitude  (km ) 

Figure  16.  Molecular  effects  variation  on  the  required  number  of  integrated  pulses. 
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Required  Number  of  Integrated  Pulses  Required  Number  of  Integrated  Pulses 


Altitude  (km) 

17.  Viewing  angle  variation  effects  on  the  required  number  of  integrated  pulses. 


Altitude  (km) 


Figure  18.  Aircraft  altitude  variation  effects  on  the  required  number  of  integrated  pulses. 
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Altitude  (km ) 


Figure  19.  Laser  line  variation  effects  on  the  required  nximber  of  integrated  pulses. 


Altitude  (km) 

Figure  20.  Required  number  of  integrated  pulses  to  maintain  a  desired  total  SNR. 
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4.3.  Range  and  Velocity  Resolution  Considerations 

Range  and  velocity  resolutions  for  a  pulsed  Doppler  lidar  system  are  inversely  related.  If 
improved  velocity  resolution  (and  accuracy)  is  desired,  then  there  is  a  sacrifice  in  the  range 
resolution  (and  accuracy). 

The  maximum  range  resolution,  AR,  is 


AR  = 


(8) 


where  tp  is  the  laser  pulse  length.  The  velocity  resolution,  AV,  is 


XAf  X 

2  -2tp’ 


(9) 


where  Af  is  the  I.F.  (intermediate  frequency)  Doppler  shift  due  to  the  change  in  line-of-sight 
velocity. 


A/=^^.  (10) 

This  leads  to  the  relationship  between  the  velocity  and  range  resolution;  that  their  product  is  a 
constant  dependent  only  upon  the  lidar  wavelength, 


ARAV  = 


cX 


(11) 


For  the  present  Ballistic  Winds  coherent  lidar,  the  constant  product  (AR  AV)  has  a  value 
of  151  m^/s.  With  a  BW  laser  pulse  length  of  0.25  ps,  the  maximum  range  resolution  is  37.5m 
and  the  velocity  resolution  is  4  m/s. 

4.4.  Velocity  Accuracy  Considerations 

The  velocity  accuracy,  Cy,  which  is  related  to  the  Doppler  frequency  accuracy  (a y),  is 
dependent  on  the  SNRy^, 


^  Skolnik,  M.I.,  editor  (1970)  Radar  Handbook,  McGraw-Hill,  New  York. 
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(12) 


A, 

=  2 

l3  1  X 
~n  2%  tpSNRv’ 

Using  tp  for  the  BW  lidar,  the  relationship  in  Eq.  (12)  is  illustrated  in  Figure  21.  The 
single  hit  velocity  accuracy  for  the  BW  simulation  is  shown  in  Figure  22  for  different  viewing 
angles.  Once  again,  molecular  backscatter  effects  are  not  included  in  the  return  signal.  The 
curve  in  this  figure  is  just  the  mirror  image  of  the  single  pulse  SNR  curves  since  the  accuracy  for 
fixed  pulse  length  depends  only  on  the  inverse  of  the  SNR. 

When  pulse  summing  is  used  to  increase  SNR,  the  effective  velocity  accuracy  reaches  a 
plateau  that  depends  on  the  required  integrated  SNR.  There  are  negligible  variations  in  the 
integrated  accuracy  from  the  various  lidar  parameter  variations  analyzed.  Figure  23  shows  the 
integrated  pulse  velocity  accuracy  with  viewing  angle  variations  analyzed.  The  velocity 
accuracy  reaches  about  0.5  m/s  (for  10  dB  power)  by  9  km  altitude. 

To  further  establish  the  validity  of  the  statement  that  the  velocity  accuracy  only  improves 
as  for  the  Ballistic  Winds  coherent  lidar,  the  analysis  of  Dunn*  was  utilized.  Dunn  gives  the 
lower  bounds  on  velocity  estimation  using  single  hit  accuracy/resolution  with  polynomial 
smoothing  and  centered  data  intervals  as. 


'  AV^' 

-1 

'l2AR^\^' 

-1“ 

= 

N 

+ 

1 

1 

1 _ 

(13) 


where  v  is  the  laser  PRF.  For  the  BW  lidar  system  values  give,  the  improvement  in  velocity 
accuracy  varies  as  for  the  BW  WSMR  lidar  system,  as  expected  for  incoherent  pulse 
summation.  The  first  term  in  Eq.  (13)  dominates  since  (M  v)^  is  much  larger  than  (AV)  .  This 
relation  is  plotted  in  Figure  24. 


*  Dunn,  K.-P.  (1987)  Lower  Bounds  on  Acceleration  Estimation  Accuracy,  MIT  Lincoln 
Laboratory,  Technical  Report  #790,  ESD-TR-87-080,  5  October,  ADA18838314. 
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Velocity  Accuracy  Improvement  Factor 


Velocity  Accuracy  (m/s) 


5.  WIND  EFFECTS  ON  AIRCRAFT 


This  section  describes  the  effects  of  the  crosswind  on  the  operation  of  the  Doppler  lidar 
and  how  the  results  of  the  tests  must  take  this  into  account.  During  the  final  leg  of  the  racetrack 
pattern  flown  by  the  B-52,  the  aircraft  is  flying  orthogonal  to  the  prevailing  WSMR  westerly 
winds.  In  order  to  fly  due  north  during  the  bomb  run,  the  pilot  must  fly  slightly  into  the 
crosswind  (F^v)  to  maintain  his  course  as  shown  in  Figure  25.  This  heading  angle  into  the  wind 
(0J,),  depends  on  the  aircraft  velocity  (F^/c)  as, 


(14) 


where  the  effective  A/C  ground  velocity  (Fg^)  is, 


Veff  =  VaIC  COSO;, 

-  V^^/c  ^w- 

This  heading  angle  correction  for  the  final  leg  is  shown  in  Figure  26  for  the  two  expected 
extremes  of  the  aircraft’s  velocity.  The  heading  angle  correction  for  the  prevailing  crosswind  is 
between  10°  and  15°. 


Figure  25.  Aircraft  orientation  during  bomb  drop  test  leg. 
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Wind  Velocity  (m/s) 


The  heading  into  the  wind  affects  the  Doppler  lidar  wind  measurements.  Given  that  the 
lidar’s  line-of-sight  is  always  orthogonal  to  the  A/C  velocity  vector,  the  LOS  wind  velocity 
^  measured  by  the  lidar  is, 

Vwh,os 

=  Vw  cosG/,  sincp;^. 


Eq.  (16)  includes  the  effects  of  the  off  nadir  viewing  angle  (914)  as  illustrated  in  Figure  27. 


Figure  27.  Off  nadir  viewing  angle  effect  on  LOS  wind  vector. 
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The  relative  LOS  velocity  reduction  with  no  nadir  effects  included  (cp^  =  90°)  is 
illustrated  in  Figure  29.  This  loss  due  to  the  A/C  heading  into  the  wind  to  maintain  course,  is 
less  than  a  4%  effect  and  is  not  a  cause  for  concern. 

When  the  off  nadir  viewing  angle  is  included,  the  normalized  LOS  horizontal  wind 
component  measured  is  further  reduced  by  the  sine  of  the  off  nadir  angle.  This  further  reduction 
is  shown  in  Figure  29.  A  15°  off  nadir  angle  reduces  the  measured  horizontal  wind  component 
to  26%  of  its  value.  Depending  on  the  magnitude  of  the  prevailing  winds,  this  can  reduce  the 
velocity  measurement  to  within  the  velocity  resolution  limit,  resulting  in  incorrectly  converting 
to  the  horizontal  wind  speed.  Increasing  the  off  nadir  angle  to  30°  almost  doubles  the  required 
components  measured  value. 


Wind  Velocity  (m/s) 


Figure  28.  LOS  velocity  reduction  due  to  crosswind. 
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Figure  29.  Relative  LOS  horizontal  wind  component  due  to  lidar  off  nadir  angle. 


6.  SUMMARY  AND  CONCLUSIONS 

This  study  gives  a  methodology  for  predicting  the  performance  of  coherent  lidar  systems 
which  use  the  atmospheric  background  aerosols  as  a  backscattering  medium  to  obtain  a 
measurement  signal.  A  real  time,  high  altitude  wind  profile  can  be  obtained  and  used  by  onboard 
aircraft  targeting  systems  utilizing  coherent  Doppler  lidars  to  make  more  refined  ballistic  wind 
corrections. 

The  major  conclusion  of  this  preliminary  study  on  the  current  Ballistic  Winds  Doppler 
lidar  is  that  the  system,  with  pulse  summing  techniques,  is  feasible  and  can  make  the  necessary 
measurements  during  aircraft  flights  at  White  Sands  or  elsewhere.  Although  the  present  analysis 
showed  that  summing  less  than  100  pulses  can  measure  the  high  altitude  winds,  there  is  one 
major  caveat.  That  is,  this  study  used  the  “nominal”  atmospheric  conditions  and  constituents 
contained  in  the  LOWTRAN7  model,  and  no  real  atmosphere  is  ever  “nominal”.  Experience  at 
the  Utah  Test  and  Training  Range^  has  shown  that  there  are  days  when  there  are  not  enough 
aerosols  present  to  obtain  measurements  back  from  even  a  powerful  groimd  based  CO2  coherent 
lidar  let  alone  a  small  aircraft-based  lidar  system.  Experimental  results  have  shown  that  the 
LOWTRAN  aerosol  models’  loadings,  at  least  at  high  altitude,  give  larger,  more  optimistic 


^  Beland,  R.R.  (1995)  Private  Communication,  Air  Force  Phillips  Laboratory/Geophysics 
Directorate  (PL/GPOA),  Hanscom  AFB,  MA. 
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signal  returns.  More  pulse  summing  will  probably  be  necessary  in  order  to  obtain  reliable  signal 
levels  from  actual  measurements  programs. 

There  are  many  methods  available  to  increase  the  signal  return  into  the  detector.  The 
most  obvious  is  to  increase  the  transmitted  pulse  energy,  but  this  is  not  always  possible  since 
most  new  applications  are  pushing  the  state-of-the-art  and  increased  energy  is  difficult  to  come 
by.  Other  improvement  techniques  are  to  increase  the  number  of  pulses  summed,  improve  the 
optics  efficiency,  operate  at  a  different  wavelength  with  more  backscattering  and  less  attenuation, 
obtain  a  more  efficient  detector/optics  coupling,  increase  the  off  nadir  look  angle,  etc.  Each  of 
these  applications  has  problems  associated  with  its  implementation. 

One  of  the  easiest  methods  to  increase  the  single  pulse  SNR  and  avoid  the  difficulty  of 
pulse  summing  is  to  increase  the  size  of  the  collection  aperture.  There  is  a  one-to-one  increase  in 
improvement  with  the  increase  in  receiver  diameter.  This  is  shown  in  Figure  30  where  the  7  cm 
collector  is  assumed  to  be  the  baseline  design.  Increasing  the  receiver  aperture  from  7  cm  to  10 
cm,  the  required  pulse  integration  factor  would  decrease  by  1.4. 


Figure  30.  SNR  improvement  factor  vs.  collection  aperture  size. 


Although  the  system  will  be  initially  installed  on  a  B-52  at  an  off  nadir  viewing  angle  of 
15°  which  is  the  maximum  angle  due  to  optics  and  space  limitations,  this  angle  is  far  from 
optimum  for  wind  measurements.  The  off  nadir  viewing  angle  should  be  as  large  as  possible, 
consistent  with  meeting  the  requirements  for  measuring  the  wind  at  the  desired  minimum 


28 


altitude.  An  angle  of  at  least  30°  would  increase  the  horizontal  wind  component  measurement  by 
a  factor  of  two. 

Molecular  absorption  is  important,  but  maintaining  the  laser  wavelength  tuned  near  the 
peak  of  a  transmission  line  is  relatively  easy.  Furthermore,  most  absorption  lines  are  broad 
enough  to  have  a  relatively  small  effect  on  the  laser’s  transmission. 

It  cannot  be  emphasized  enough  that  these  analyses  were  made  with  LOWTRAN7 
standard  mid-latitude  spring/summer  atmospheric  profiles.  The  real  atmosphere  can  be  (and 
usually  is)  appreciably  different  -  orders  of  magnitude  variation  in  attenuation  and  scatter 
coefficients  compared  to  “nominal”  are  not  uncommon.  This  leads  to  the  open  question  of  what 
the  aerosol  concentration  at  the  high  altitudes,  especially  the  tropopause,  really  is  at  the  sites  of 
interest.  The  tropopause  is  a  highly  variable  region  and  can  greatly  affect  the  SNR  and  pulse 
summing  predictions.  For  any  reliable  predictions  of  field  test  results,  actual  atmospheric 
profiles  especially  those  relating  directly  to  the  laser  backscatter  coefficient  should  be  measured 
at  the  site. 
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